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NaTIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
ADVANCE CONFIDENTIAL REPORT 

FLIGHT INVESTIGATION OF TEE VARIATION OF DRAG COEFFICIENT 
WITH MACH NUMBER FOR THE BELL P-39N-1 AIRPLANE 
By Welko E. Gasich and Lawrence A„ Clousing 

SUMMARY 


An investigation of the effect of compressibility on the 
drag of a Bell P-39N-1 airplane has been made in flight as 
part of a drag study on several high-speed airplanes currently 
in production. The Nach number range covered in these tests 
was from 0.2 to about 0,8. 

The minimum drag coefficient at low Mach number was found 
to be 0.022. The drag coefficient started to increase at a 
Mach number of 0,62, reaching a value almost three times as 
great at a Mach number of 0.80. From the results of the tests 
it appears that the terminal Mach number for the airplane dur- 
ing a dive from service ceiling is about 0,30, 


INTRODUCTION 


One of the compressibility effects which occurs on air- 
planes in flight at speeds approaching those of the speed of 
sound is that of a large increase in the airplane drag coef- 
ficient. The Mach number at which this drag increase occurs, 
and the ability to predict its magnitude, is a subject of in- 
tense research at present, because the high-speed performance 
of modern airplanes is affected by the nature of the drag in- 
crease. Also, the terminal dive speed of airplanes is deter- 
mined to a large extent by the nature of the increase of drag 
at high Mach numbers. At present, reliable full-scale data 
on the drag characteristics of airplanes at high Mach numbers 
are quite limited. 
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As a result, the Ames Aeronautical Laboratory has under- 
taken drag- investigations in flight of several high-speed air- 
planes currently in production* The present investigation was 
conducted to determine the manner in which the drag coeffi- 
cient of a Bell P-39N-1 airplane varied with Mach number* 
Measurements were taken up to the highest value of Mach number 
attainable* 

In conducting the investigation, it was found necessary 
to devote special study to the technique for the determination 
of drag characteristics of airplanes in flight. A description 
is presented herein of the methods used, and the results ob- 
tained through the use of different methods of computing the 
drag are compared. 


SYMBOLS 


Tn the derivation of the formulas by which the drag coef- 
ficient is obtained, the following symbols are used: 

F total force parallel to flight path acting to accelerate 
the airplane, pounds 

s distance through which force acts, feet 
m mass of airplane, I7/g, pound-seconds squared p>er foot 
W weight of airplane, pounds 

a acceleration of airplane along flight path (in direction 
of F) , feot per second squared 

g acceleration due to gravity, feet per second squared 

h true altitude, feet 

V true airspeed, miles per hour 

D total drag of airplane, pounds 

q dynamic pressure, pounds per square foot 

S wing area, square feet 
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T total thrust, pounds 

the algebraic sum of the components along the airplane 
X-axis of the airplane acceleration and the accelera- 
tion due to gravity in terms of the standard gravi- 
tational unit (32.2 ft/sec' J ), positive v/hen directed 
forward 

Ax longitudinal acceleration factor due to attitude of air- 
1 plane in terrestrial gravitational field 

Ax o longitudinal acceleration factor due to acceleration of 
airplane through space 

Ag the algebraic sum of the components, along the airplane 
Z-axis, of the airplane acceleration and the accel- 
eration due to gravity, in terms of the standard 
gravitational unit (32.2 ft/sec 2 ), positive when di- 
rected upward 

normal acceleration factor due to attitude of airplane 
in terrestrial gravitational field 

A^ normal acceleration factor due to acceleration of air- 
s plane through space 

8 angle of flight path as measured from horizontal, degrees 

CX angle of attack of thrust line, degrees 

M Mach angle 


DESCRIPTION OF THE AIRPLANE 


The Bell P-39N-1 airplane is a single-place, low-wing, 
cantilever monoplane powered by a 1200 brake horsepower (take- 
off rating) Allison V-17 10-85 liquid-cooled engine driving a 
three-blade Aeroproducts propeller. Figure 1 is a general- 
arrangement drawing of the airplane. Figure 2 is a photograph 
of the airplane as instrumented for the flight tests. The 
following specifications of the airplane were taken almost en- 
tirely from references 1 and 2: 
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Airplane, general 

Span . . . . 

Length 

Weight (normal end approx, as flown) 

Center of gravity (for normal gross 
weight and approx, as flown). . . 


34.0 ft 
30.167 ft 
7629 lb 

0.285 M.A.C. 


Wing 

Airfoil section, root ........ NASA 0015 

Airfoil section, tip NACA 23009 


Area ... ...... •••••»•• 213.2 sq ft 

Engine 

Type ................. Allison V-1710-85 

Ratings ( bhp/ rpm/ altitude ) 

Take-off .....o. .«»»•« 1200/3000/s ea level 

Military ............. 1125/3000/15,500 

Normal 1000/2600/l4-,000 


Gear ratio .............. 2.23:1 

Propeller 


Type 


Aeroproducts, 
const ant- speed, 
ho llow- shaft 


Blade design ..... 

Number of blades ..... 

Activity factor/blade ........ 

Thickness ratio, 75-percent radius . . 
Diameter ...«.o.... .o. • * 


A-20-156-17 


3 

98.5 

0.08 

11 ft, 7 in. 
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INSTRUMENTATION 


Standard NACA instruments were used tc record photograph- 
ically ^ as a faction of time, quantities from which the fol- 
lowing variables could be obtained: indicated airspeed, pres- 

sure altitude, normal acceleration, longitudinal acceleration* 
engine manifold pressure* engine speed, and approximate angle 
of attack of the thrust line* The instrument for measuring 
longitudinal acceleration was made especially sensitive for 
the purpose of the tests* 

The recording instruments, as installed in the airplane, 
could be read to £2 miles per hour for the indicated airspeed* 
± 200 feet for the altitude, ±0 o lg for the normal acceleration* 
±0 o 01g for the longitudinal acceleration, ±0*2 inch of mercury 
for the manifold pressure, ±20 rpm for the engine speed, and 
±0^2° for the angle of attack* 

The temperature was measured by the pilot f s service ther- 
mometer o Temperature surveys were made in ascending and de- 
scending flight in an effort to eliminate errors in tempera- 
ture reading caused by the lag effect in the system* The 
average temperature read was corrected for the temperature 
rise caused by the compression of the air due to the spcod of 
the airplane* 

The head used for the airspeed and altitude measurements 
was freely swiveling and was mounted on the end of a boom ex- 
tending about 4 feet ahead of the leading edge of the right 
wing at a spanwisc location about 7 feet inboard of the right 
wing tip 0 The airspeed head consisted cf two separate static- 
pressure tubes (one of which was connected to the airspeed re- 
corder and the other to the altitude recorder) with a single 
total-pressure tube located between them-, The airspeed and 
altitude recorders were mounted in the right wing at the base 
of the booiiie The total- and static-pressure linos to the air- 
speed recorder were balanced so that a sudden pressure change 
equally apolied to both total- and static-pressure linos at 
the airspeed head caused no reading of the airspeed recorder v 
The lines to the airspeed recorder and altitude recorder wero 
of about the same length and size, being about 5 feet long 
and of 0*12 inch inside diameter* The magnitude of the lag 
effects in the altitude and airspeed recording systems was 
measured by a specially built ground setup <s ..and v/as found in 
both systems to be smaller than an error in pressure equiva- 
lent to 5 feet of altitude with the airplane in a dive at 
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terminal Mach number* The recording and service static heads 
were calibrated for position error by comparing the readings 
of the recording and service altimeters with the known pres- 
sure altitude as the airplane was flown at several speeds past 
a reference heights It was assumed that the correct total 
pressure was obtained a Calibration of the recording head in 
the Ames 16-foot wind tunnel showed that the error in recorded 
airspeed, due to the difference in Mach number between the 
highest value obtained in the flight calibration (0 o 50) and 
the highest value obtained in the flight tests (0»80), was less 
than 1 percent. Indicated airspeed, as used in this report, 
was computed according to the formula by which standard air- 
speed meters are graduated (gives true airspeed at standard 
sea-level conditions). The formula may be written as follows: 


V i = 1703 




1/2 


where 


V correct indicated airspeed, miles per hour 

H free-stream total pressure 

p free-stream static pressure 

p n standard atmospheric pressure at sea level 

The angle of attack of the thrust line was indicated by 
a vane mounted on the forward end of a boom located at a span- 
wise station 6.7 feet inboard from the left wing tip an d ex- 
tending 3.3 feet ahead of the wing leading edge similar to the 
airspeed-head installation. The angle of attack as measured 
by th* s vane was used in this analysis without any corrections 
being applied for position error. 


METHODS OF ANALYSIS 


It was found necessary in the course of these analyses, 
as will bo explained later, to calculate drag coefficient from 
high-speed-dive data for Mach numbers above 0.5 and from speed- 
power data for Mach numoers below this figure. Three differ- 
ent methods wore applied to the dive date, while a single meth- 
od was applied to the speed-power data. 
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The methods ■which were applied to the dive data will be 
referred to in this report as (1) force method, (2) energy 
method, and (3) acceleration method* The equations derived in 
appendix A are as follows: 

(l) Force method 


T cos a - W 


Cp e 


dh/dt dV/dt 


- V 


qS 


_I> r 


(2) Energy method 


T cos a - ^ — ' h + ~ 

C D » — X 2S-: 

qS 


(3) Acceleration method 


CD 


T cos a + W(Az sin a - A^ cos a ) 

~qs~ 


The equations shown in the foregoing for the force and energy 
methods are identical except for form 0 It was felt that a dis- 
tinction should be made, however, in order to emphasize the 
significance of the terms within the brackets* Also, the im- 
plied differentiations were advered to in applying the equa- 
tions* 


The foregoing methods yielded poor results when applied 
to test data where the Mach number was below 0*5; therefore, 
another method was developed which was applicable to data ob- 
tained during speed-power runs in level flight© This method 
is discussed in appendix B, but is not completely derived be- 
cause of the rather lengthy series of computations involved. 

In evaluating the total thrust developed by the engine- 
propeller combination, it was necessary to include both the 
exhaust thrust and propeller thrust. The exhaust thrust was 
estimated by the method used in reference 3 and the propeller 
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thrust was determined from estimations of propeller efficiency 
and engine power. Propeller efficiency was estimated by the 
method described in appendix C using charts and data presented 
in references 4, 5, and 6. 

The engine brake horsepower was determined from the manu- 
factxircr's engine-power charts by entering readings of the re- 
cording altimeter, rpm recorder, and manifold pressure recorder. 
In many instances, the manifold pressure at given altitude and 
engine- speed conditions exceeded that shown by the engine-power 
charts. In such cases it was assumed that a line of constant 
rpm and manifold pressure could be extrapolated linearly to the 
higher altitude to determine engine power# It might be pointed 
out that this assumption is not exactly true, because the ad- 
ditional manifold pressure is not supplied by the gear-driven 
supercharger within the engine but by the ram effect in the fuel- 
induction system; hence, the manifold pressure is not obtained 
entirely by the expenditure of power in the supercharger,. The 
amount of ram pressure is, however, such a small percentage of 
the total manifold pressure that it was believed that no seri- 
ous errors resulted from this assumption. 


TESTS, RESULTS, AND DISCUSSION 


The data were obtained in high-speed dives and in speed- 
power runs with the airplane in the clean condition with oil 
and coolant shutters one-half open (flush with the external 
lines of fuselage). Pressure orifices wore present in the 
empennage and right wing , having been placed there for other 
test purposes, and an extra aerial mast was mounted ahead of 
the cockpit canopy. A research angle-of- attack head and boom 
was mounted on the left wing and a research airspeed head and 
boom was mounted on the right wing. A bomb rack and sway 
braces were mounted under the center of the fuselage. In or- 
der to insure that the landing- wheel doors did not open at 
high speeds and pull-outs, special latches were installed so 
that the doors wore positively locked during flight, A list 
of the dives involved in these tests and the approximate 
flight conditions during the dives are presented in table I. 

In order that the force, energy, and acceleration methods 
of determining the drag coefficient could be compared in re- 
gard to accuracy and consistency, two dives, the time his- 
tories of which are shown in figures 3(a) and 3(b), were ana- 
lysed by each of the methods# 
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Figure 3(a) shows a time history of a power-on dive (dive 
No. 1 of table i) from 30*000 feet at full throttle and 3050 
rpm. The maximum indicated airspeed reached was 464 miles per 
hour at 12,500 feet. The maximum Mach number reached in the 
dive was 0.765, The airplane was pulled out at 11,500 feet 
at a normal acceleration of 6.0g. It should be noticed that 
the manifold pressure kept rising constantly until a value of 
50 inches of mercury was reached, at which point there was a 
sharp break and the pressure became constant. This break was 
due to the pilot f s throttling back in order to keep the mani- 
fold pressure below the maximum permissible value of 51 inches 
of mercury. 

Figure 3(b) shows time history of a power-on dive (dive 
No, 2 of table i) from 28,000 feet at part throttle and 2800 
rpm. The maximum indicated airspeed reached in this dive was 
471 miles per hour at 10,300 feet. The maximum Mach number 
reached was 0.777. The reason for the sudden decrease in manifold 
pressure at the beginning of the dive is due to loss of power 
caused by malfunctioning of the fuel-induction system of the 
engine when negative accelerations are incurred. 

The airplane drag coefficient calculated from the data 
obtained in the dive shown in figure 3(a)* as determined by 
the three methods of evaluation, is shown in figure 4(a) with 
airplane drag coefficient Cq plotted as a function of Mach 
number M, The results are presented only for Mach numbers 
above 0,6, The variation of Cq with M as determined by 
the acceleration method appears to be the most reasonable 
variation of those shown. The energy and force methods are 
consistent within themselves to a small degree, but are in no 
way consistent with the acceleration method. The drag-coef- 
ficient curves determined by the force and energy methods, 
in fact, show an illogical variation with Mach number. 

An attempt was made to investigate the reason for this 
variation by assuming that the curve of Cp plotted as a 
function of Mach number as determined by the acceleration 
method was correct. From this curve the data were worked 
backward until a time history of the true airspeed was deter- 
mined. A comparison of the actual variation of airspeed with 
time that was obtained from the data showed that in no case 
was there a difference of more than 5 miles per hour true 
airspeed. This error in measurement of true airspeed is en- 
tirely possible, especially when it is realized that in deal- 
ing with true airspeed an accurate temperature measurement is 
necessary to obtain true airspeed from indicated airspeed. 
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Hence, it would appear that inherently the energy and force 
methods are too sensitive to minor errors in the measurement 
of airspeed to make the methods usable in flight work. The 
inconsistency between force and energy methods may be due to 
the fact that for the force method two graphical differenti- 
ations are necessary, one to determine the dive angle and the 
other to determine the longitudinal acceleration, while only 
one differentiation is necessary for the energy method,. 

The drag results of the dive of figure 3(b) are shown 
on figure 4(b) „ In this case the results of the energy and 
force methods do not agree with each other as well as in the 
previous example, although in this case the energy and force 
methods give lower drag coefficients at Mach numbers less 
than 0.71 as compared to values presented in figure 4(a). 

Because of the apparent inconsistency of results obtained 
by the energy and force methods, all the dive data we re finally 
derived by the acceleration method,, These data are shown in 
figure 5. The results are consistent and the test points de- 
termine a well-defined curve of drag coefficient as a function 
of Mach number for a lift-coefficient range of 0.03 to 0.09 « 

The greatest inaccuracy in computed drag coefficient at the 
high Mach numbers should not exceed 2 percent, because the en- 
gine thrust (propeller and exhaust) does not exceed 7 percent 
of the total drag at these Mach numbers. The possible error 
rises, however, at low Mach numbers, due to the value of the 
thrust more closely approaching that of the drag. At a Mach 
number of 0.5, the thrust is equal to the drag. 

The drag coefficients at Mach numbers below 0.5 have been 
calculated by the method outlined in appendix B and are shown 
in figure 5, To obtain data for this method, the airplane was 
flown in straight level flight at various power conditions at 
an altitude of 15,000 foet, thus obtaining a speed-power curve 
for the airplane. The results of the data after they were re- 
duced and plotted as figure 6 indicate an airplane drag co- 
efficient of 0.022 at a lift coefficient of zero, which seems 
to be a reasonable figure, since the particular airplane tested 
was aerodynomically rather unclean. 

It is of interest to note in figure 5 that one of the 
points plotted was obtained from dive data during a time in 
the dive recovery in which the Mach number was decreasing after 
having reached a high value. The position of this point above 
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the general curve may possibly be attributed to a hysteresis 
ef* i ect caused by the continuance of compressibility stall on 
some portion of the airplane following flight at high Mach 
number,, It may also have been caused by increased bulging 
or distortion of the surface of the airplane as it dived 
into the lower atmosphere where a larger value of dynamic 
pressure occurs at a given Mach number than is the case at 
high altitude© 

Data are presented in figure 5 on the section critical 
Mach numbers (the flight Mach number at which the local air- 
speed over the surface reaches sonic velocity) for various 
spanwise stations on the wing of the airplane* The various 
spanwise stations at which the section critical Mach number 
of the wing was determined are shown in figure 7 0 The varia- 
tion of wing thickness with span is also shown in this figure* 
The values of critical Mach number of each of the spanwise 
stations shown in figure 7 were determined as described in ap- 
pendix D and are shown in figure 8© The pressure-distribution 
measurements used in determining the critical Mach number of 
each spanwise station were obtained during the flight tests 
reported in reference 7, and are on file at Ames laboratory. 

The critical Mach number for the four approximately symmetri- 
cal inboard sections (0015 to 0012) is about 0 o 68© The crit- 
ical Macn number for the tip section (23010) is much lower 
than that for the symmetrical sections, because the tip sec- 
tion is operating at a lift coefficient less than the optimum 
lift coefficient for maximum critical Mach number e The values 
of wing critical Mach number are spotted on figure 5 to show 
their relation to the airplane Mach number of drag divergence. 
(Drag divergence is defined as the point at which the drag coef- 
ficient begins to increase over its low speed value.) 

It is seen that the drag coefficient starts increasing at 
a Mach number of 0,62, about 0.03 higher than the tip-section 
critical Mach number. The critical Mach number on the rest of 
the wing is not reached, however, until approximately 0.06 Mach 
number after the airplane Mach number of drag divergence is 
reached* It is probable that shock waves developing on parts 
of the airplane (e.g«, canopy, duct entrances, etc.), other 
than the wing tip, contribute to the early drag increase. 

It is interesting to notice that the value of dC^/dM keeps 
gradually increasing up to a Mach number of 0.75. The value of 
dCp/dM between Mach numbers of 0.75 and 0.80 is approximately 

0.5. 
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The highest Mach number reached in any dive was about 0*8 
which appears to be the terminal Mach number for the Bell P-39N-1 
airplane when dived power on from its service ceilings No dif- 
ference in the maximum Mach number attainable was discerned be- 
tween dives where the entry was made by rolling over into the 
dive and those in which entry was made by pushing down into the 
dive. Calculations of the terminal Mach number which could be 
reached in a vertical dive from 34,000 feet to sea level were 
made by a step-by-step process using the drag curve of figure 5 
with a slight extrapolation. Since at high Mach numbers the pro- 
peller thrust in powor-on flight is avery small proportion of 
the total thrust component, the effect of the propeller thrust 
was neglected in these computations. The results indicated that 
the highest Mach number would be obtained at an altitude of about 
19,000 feet and that this Mach number would be 0*802, which is 
but slightly higher them the maximum Mach number obtained in the 
drag tests.. The computations indicated a decrease in Mach number 
as the airplane continued the dive below 19,000 foot* 

The highest Mach number reached in level flight during the 
tests wa^ G c 4£u Data of reference 8 show, however, that with 
war emergency power a Mach number of 0*54 can be reached in level 
flight with a Bell .P-3917-1 airplane* With a Mach number of 
drag divergence of 0*62, the airplane could, therefore, be 
flown about 60 miles per hour faster than at present oefore 
the compressibility effect on drag would start to limit the high 
speedy 


CONCLUSIONS 


1 0 The minimum drag coefficient at low Mach number for 
the Bell P-39N-1 airplane was found to be approximately £).022 o 

2, The Mach number of drag divergence (that Mach number at 
which the drag coefficient started to increase from its low- 
speed value) was 0*62 0 

3, The r.rximum Mach number attained in the course of the 
tests was about 0.80, which appears to be the terminal Mach 
number of the airplane. At this Mach number, the drag coeffi- 
cient was about 0.060. 

4, The slope of the curve of drag coefficient as a func- 
tion of Mach number (dC^/dM) in the high Mach number region 
was about 0.5, 
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5* The critical Mach numbers of the different airfoil sec- 
tions of the wing as determined from pressure-distribution 
measurements are, in general, in good agreement with the theo- 
retical values of critical Mach number for the various airfoil 
sections. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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APPENDIX A 


Three methods are derived by which the airplane drag coef- 
ficient Cp may be obtained from flight-test data# The three 
methods are referred to as (l) the force method, (2) the energy 
method, and (3) the acceleration method# 


Force Method 

The drag of the airplane may be determined by use of rec- 
ords of altitude and airspeed,, Equating the forces along the 
flight path (fig# 9) gives 

F - T cos o& + W sin 0 - D 


or 


from which 


— a = T cos a+W sin 0 — D 
g 


D * T cos a + W sin 0. - 

a = d7/dt 

sin 0 = ~ 

7 



Substituting for sin 0 and 


a 


or 



where sin 0 and a are determined by graphical differentiation 
with respect to time, of altitude and velocity along the flight 
path# 
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Energy Method 

This method involves the equating of the changes in po- 
tential and kinetic energy to the product of force and dis- 
tance as follows: 


since 


and 


then 


and 


d(work) z Fds * d(mgh) + d --- 

E/ 

= W dh + W 


- W dh + 


V dV\ 

* ) 


ds “ V dt 


F « T cos a - D 


r 


(T cos a - D) V dt ■ W ^ dh + V ^ 

D s T cos a - I A ( h + 

V dt \ 2g/ 


Cd = 


T cos a - •& 

V dt 

as 


t£> 


It is seen that this method does not require the use of 
an accelerometer. The only instruments required are an al- 
timeter and an airspeed recorder. The graphical differentia- 
tion of the quantity (h + V /2g) with respect to time is 
required to evaluate the latter part of the right-hand term 
of the above equation. This equation is essentially the same 
result as is obtained by the force method, except that one 
graphical differentiation is required instead of two. The re- 
sults obtained from the two methods differ only by errors due 
to the work-up of the basic flight data. 
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Acceleration Method 

Determination cf the drag of an airplane in flight may 
bo made by consideration of the accelerations involved e 

As before. 


D a T cos a + W(sin 0 - a/g) 

sin 9 is measured by the components of the accelerometer due 
to their orientation in the gravitational fieldo Therefore 


sin 0 » Ay cos a.+ A? sin a 

1 1 

a/g is measured by the components of the accelerometer due 
to their acceleration along the. flight patho Therefore 


a 

i 


Ay cos a 

A 3 


A 7 sin a 
L 3 


hence 

D s T cos a+ W(A X cos a - A x cos a + A^ sin a - A^ sin a) 
J ^ 2 1 2 

and since 


A X a ' A X, 


A z “ A z x - A z a 


then 


or 


D a T cos a + W(Ag sin a 


CD 


T cos a + ^(A? sin a 

qS 


- A x cos ct ) 

- A x cos a) 
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APPENDIX B 


The following discussion presents the method by which the 
drag coefficient at low Mach numbers was obtained from speed- 
power data. The method is based on the assumption that the 
airplane polar may be represented by a parabola in the normal 
flying rangej for example, lift coefficients from 0„0 to 0,8, 

The equation for the parabola is then written as 


where 


CD = c Dp + 


cl 3 

v eA 


Cp total airplane drag coefficient 

Cp^ effective parasite drag coefficient 

e airplane efficiency factor 

A aspect ratio 

By proper mathematical manipulation, the following equa- 
tion for indicated thrust horsepower required for level flight 
may be determined: 


th Pi 


thp(a) 1 ,z = ' 6,82 


lif f + 0*332 ( i] 

100/ eVi V b / 


where 


f 

C n S, sq ft 


Vi 

^true, mph 


w 

airplane gross 

weight 

b 

airplane span. 

ft 


The above-mentioned formula for the power-required term de- 
fines a single curve of indicated power required as a function 
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of indicated airspeed for all altitudes. By the proper appli- 
cation of the foregoing formula, flight-test data may be re- 
duced to give the airplane parameters e and f. This may be 
accomplished by representing the power-required equation as 
linear; that is, represented by the intercept equation of a 
straight line 


where 


and 



y = e 

Hence the power-required equation becomes 


6,82 (Vi/100) 3 f 0,332 (w/b) 2 

ammm — ~ m. mm ~ mm —mm mm mm mmmm mi mm Mi mm mm mm h mi 

thpi eV i thpi 


If the values of Sj*82—iYi/i22) are plotted against values of 
3 th Pi 

2 then the intercepts of the line passing through 

Vi thpi 


the given test points determine the equivalent flat-plate area 
f and airplane efficiency' factor e« In order to determine 
the indicated thrust power, it is necessary to define the en- 
gine brake horsepower by engine charts or a torquemeter, as the 
case may be c , From the power conditions (rpm, altitude, etc*) 
the propeller efficiency may bo determined and hence the pro- 
peller thrust power* To this power must be added the exhaust 
thrust power v if any 0 Exhaust thrust may be estimated by the 
method of reference 3 0 By multiplying this value of thrust 
power by Jo, the indicated power required (thp-j_) is de- 


termined, leaving only the gross weight and indicated airspeed 
to be evaluated * 
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APPENDIX C 


Propulsive efficiency was estimated from charts presented 
in reference 4 corrected for tip compressibility effects using 
the data obtained from reference 5 0 These data are presented 
in figure 10® The curves show the tip-speed factor as a 

function of tip Mach number for various advance-diameter 

ratios J© In using these curves the same critical tip Mach 
number was assumed to apply to the propeller on the Bell 
P-39N-1 airplane as applied to the test propeller of reference 
5, since it was of the same thickness ratio and both propellers 
had high-critical-speed tip sections* 

The compressibility losses over the blade root were esti- 
mated by the method developed in reference 6* The root cor- 
rections were then applied to the propeller efficiencies of 
reference 2 for power coefficients of 0.2 and 0*3, and for pro- 
peller tip Mach numbers of 0*8 and 1.2. These curves are shown 
in figure 11 as h r# which is the ratio of propeller efficiency 

with root losses to propeller efficiency with no root losses 0 
The over-all efficiency of the propeller was then calculated by 
multiplying the values of efficiency obtained from the chart of 
figure 10 and the chart of figure 11. 


APPENDIX D 


To establish the critical Mach number of the Bell P-39N-1 
airplane wing, use was made of normal-force-distribution data 
at five spanwise stations along the wing. These data had been 
measured in flight in connection with a separate investigation 
on the airplane. The critical Mach number was determined at 
each section for the particular section lift coefficient corre- 
sponding to an airplane lift of 0.06, which is an average lift 
coefficient for the drag curves of figure 5. Since only normal- 
f orce-distribution data were taken, the flight data had to be 
reduced to pressures on upper and lower surfaces. This was done 
as is outlined in reference 9, which presents a method for the 
rapid calculation of the pressure distribution over an airfoil 
section when the normal-force distribution and the pressure dis- 
tribution over the base profile (i.e., the profile of the same 
airfoil if the camber line were straight and if the resulting 
airfoil were at zero angle of attack) are kno*wn* Since there 
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were no pressure-distribution data available for base profiles 
corresponding to the test sections, theoretical base-profile 
pressure-coefficient distributions were used as determined from 
the above-mentioned reference. 

The resulting pressure distributions gave results which 
checked very closely the value of critical Mach number as de- 
termined by reference 10. The peak pressures at spanwise sta- 
tions A, B, C, and D occurred between 10 and 20 percent chord, 
hence the peak pressure was well defined. For station E, 
which is approximately an NACA 23010 section, the peak pres- 
sure occurs over an extremely limited portion of the airfoil 
chord near the leading edge, and the relatively few orifices 
near the leading edge prevented the peak pressures from being 
accurately established. A value for the critical Mach number 
for the tip section (station E) was arrived at, however, by 
use of reference 10. Since the pressure-distribution data 
from flight gave results that were in very close agreement with 
the theoretical data for stations A, B, C, and D, it is thought 
that the value of section critical Mach number for station E 
is probably not greatly in error 0 
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TABLE I. -TABULATIONS OF FLIGHT CONDITIONS DURING 
HIGH-SPEED DIVES OF BELL P-39N-1 AIRPLANE 


Dive 

num- 

Altitude 

! Maximum 
j normal 
acceler- 
! at ion 
| factor at 
i pull-out 

U (?) 

| Maximum 
;Mach 
| number 

1 attained 

! 

! 

Engine 

conditions 

ber 

i 

j At start 
, of dive 

i (ft) 

After 

pull-out 

(ft) 

rpm 

| Throttle 
setting 

i 

30,000 

11,500 

e.o 

0.765 

3050 

Varied 

2 

28,000 

10,300 j 

I j 

4.4 

.777 

2800 

Part 

3 

34,000 

12,100 

5.8 

j 

.797 

3000 

Full 

4 ! 

1 

32,000 

22,000 

5,0 

,715 

3000 

Part 

5 

32,000 

12,800 

7.0 

.785 

2800 

Part 

I 

6 i 

t 

30,000 

10,000 

4.6 

n 778 

2600 

Varied 


NOTE: Oil and coolant shutter position for all dives; 

one-half open (flush with external lines of fuselage). 
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Fig. 1 





Figure 1.- Three-view drawing of the Bell P-39N-1 airplane. 



Figure 2.- Three-quarter front .view of the Bell P-39N-1 airplane as instrumented 
for flight tests. 
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Figure 3 — Time history of a power-on dive from 


30.poo feet - Clean condition ^ oil and coolant shutters 
flush . B e// P~33 /VV airplane . 
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Fig. 3b 



Figure 3 , - Concluded. dell P-39M-! curp/ane. 




Ftyure 4 .*” (Comparison of | /Vie o f i u.trp/aire'dta-ft 

coeff creat with Mac h number «.y computed by. use 
©f #j 4 -force; energy end <*«< ? ehruttdn ■ m « pjie^s fVow 
?h ea5o , re»ie?7f‘a fife/] dt/zr/Viy /wfe? 
u<|rpilatie : . 
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Figure 5. - Variation of airplane dmf coefficient \knth Much 

-number as commuted ft'ow /e vef-iff/gdtdeta akd ffdiri d&fci 
faketi during six dives, Bdll P-33M-I airplane. 
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Figs,. 6,T 
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Fiyure.6.- R<^%u/is of ^p&ed- pt>^er jafo_ run af 15/jOQ.Aesi as~ piaffed -fa 
H e ter m /'?>g the dxojj ao e_f ^ Jjc.ua ?’ f o£___fjhe. & eJJ P-'Q9dVri a ne^ 



Section critical Wl&ch number 
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Figure 8- Variation of section critical Mach number with span on 
the P-39/VW airplane. 



Figure 9 r Diagram of forces umd acce/erat ions 
acting at the center of gravity. &en p-sw-i eirpt&nc. 
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figure/ /.- Chart for correofm^ propef/er efficiency 
for roof losses on Atroproducfs A-ZO-I56+-/ 7 propeVe* 
Us installed or, IZ ef ! P- 3SW-/ airplane* 
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